Independent loss events of a functional tetherin gene in galliform birds.

IMPORTANCE: Birds represent important hosts for numerous viruses, including zoonotic viruses and pathogens with the potential to cause major economic losses to the poultry industry. Viral replication and transmission can be inhibited or blocked by the action of antiviral restriction factors (RFs) encoded by the host. One well-characterized RF is tetherin, a protein that directly blocks the release of newly formed viral particles from infected cells. Here, we describe the evolutionary loss of a functional tetherin gene in two galliform birds, turkey (Meleagris gallopavo) and Mikado pheasant (Syrmaticus mikado). Moreover, we demonstrate that the structurally related protein TMCC(aT) exerts antiviral activity in several birds, albeit by a mechanism different from that of tetherin. The evolutionary scenario described here represents the first documented loss-of-tetherin cases in vertebrates.

Unraveling the palindromic and nonpalindromic motifs of retroviral integration site sequences by statistical mixture models.

A weak palindromic nucleotide motif is the hallmark of retroviral integration site alignments. Given that the majority of target sequences are not palindromic, the current model explains the symmetry by an overlap of the nonpalindromic motif present on one of the half-sites of the sequences. Here, we show that the implementation of multicomponent mixture models allows for different interpretations consistent with the existence of both palindromic and nonpalindromic submotifs in the sets of integration site sequences. We further show that the weak palindromic motifs result from freely combined site-specific submotifs restricted to only a few positions proximal to the site of integration. The submotifs are formed by either palindrome-forming nucleotide preference or nucleotide exclusion. Using the mixture models, we also identify HIV-1-favored palindromic sequences in Alu repeats serving as local hotspots for integration. The application of the novel statistical approach provides deeper insight into the selection of retroviral integration sites and may prove to be a valuable tool in the analysis of any type of DNA motifs.

Identification of GC-rich LAT genes in birds.

Linker for activation of T cells (LAT) plays a key role in T-cell antigenic signaling in mammals. Accordingly, LAT orthologues were identified in the majority of vertebrates. However, LAT orthologues were not identified in most birds. In this study, we show that LAT gene is present in genomes of multiple extant birds. It was not properly assembled previously because of its GC-rich content. LAT expression is enriched in lymphoid organs in chicken. The analysis of the coding sequences revealed a strong conservation of key signaling motifs in LAT between chicken and human. Overall, our data indicate that mammalian and avian LAT genes are functional homologues with a common role in T-cell signaling.

Cannabinoid receptor 2 evolutionary gene loss makes parrots more susceptible to neuroinflammation.

In vertebrates, cannabinoids modulate neuroimmune interactions through two cannabinoid receptors (CNRs) conservatively expressed in the brain (CNR1, syn. CB1) and in the periphery (CNR2, syn. CB2). Our comparative genomic analysis indicates several evolutionary losses in the CNR2 gene that is involved in immune regulation. Notably, we show that the CNR2 gene pseudogenized in all parrots (Psittaciformes). This CNR2 gene loss occurred because of chromosomal rearrangements. Our positive selection analysis suggests the absence of any specific molecular adaptations in parrot CNR1 that would compensate for the CNR2 loss in the modulation of the neuroimmune interactions. Using transcriptomic data from the brains of birds with experimentally induced sterile inflammation we highlight possible functional effects of such a CNR2 gene loss. We compare the expression patterns of CNR and neuroinflammatory markers in CNR2-deficient parrots (represented by the budgerigar, Melopsittacus undulatus and five other parrot species) with CNR2-intact passerines (represented by the zebra finch, Taeniopygia guttata). Unlike in passerines, stimulation with lipopolysaccharide resulted in neuroinflammation in the parrots linked with a significant upregulation of expression in proinflammatory cytokines (including interleukin 1 beta (IL1B) and 6 (IL6)) in the brain. Our results indicate the functional importance of the CNR2 gene loss for increased sensitivity to brain inflammation.

The Discovery of Chicken Foxp3 Demands Redefinition of Avian Regulatory T Cells.

Since the publication of the first chicken genome sequence, we have encountered genes playing key roles in mammalian immunology, but being seemingly absent in birds. One of those was, until recently, Foxp3, the master transcription factor of regulatory T cells in mammals. Therefore, avian regulatory T cell research is still poorly standardized. In this study we identify a chicken ortholog of Foxp3 We prove sequence homology with known mammalian and sauropsid sequences, but also reveal differences in major domains. Expression profiling shows an association of Foxp3 and CD25 expression levels in CD4+CD25+ peripheral T cells and identifies a CD4-CD25+Foxp3high subset of thymic lymphocytes that likely represents yet undescribed avian regulatory T precursor cells. We conclude that Foxp3 is existent in chickens and that it shares certain functional characteristics with its mammalian ortholog. Nevertheless, pathways for regulatory T cell development and Foxp3 function are likely to differ between mammals and birds. The identification and characterization of chicken Foxp3 will help to define avian regulatory T cells and to analyze their functional properties and thereby advance the field of avian immunology.

Dynamic Evolution of Avian RNA Virus Sensors: Repeated Loss of RIG-I and RIPLET.

Retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated protein 5 (MDA5) are key RNA virus sensors belonging to the RIG-I-like receptor (RLR) family. The activation of the RLR inflammasome leads to the establishment of antiviral state, mainly through interferon-mediated signaling. The evolutionary dynamics of RLRs has been studied mainly in mammals, where rare cases of RLR gene losses were described. By in silico screening of avian genomes, we previously described two independent disruptions of MDA5 in two bird orders. Here, we extend this analysis to approximately 150 avian genomes and report 16 independent evolutionary events of RIG-I inactivation. Interestingly, in almost all cases, these inactivations are coupled with genetic disruptions of RIPLET/RNF135, an ubiquitin ligase RIG-I regulator. Complete absence of any detectable RIG-I sequences is unique to several galliform species, including the domestic chicken (Gallus gallus). We further aimed to determine compensatory evolution of MDA5 in RIG-I-deficient species. While we were unable to show any specific global pattern of adaptive evolution in RIG-I-deficient species, in galliforms, the analyses of positive selection and surface charge distribution support the hypothesis of some compensatory evolution in MDA5 after RIG-I loss. This work highlights the dynamic nature of evolution in bird RNA virus sensors.

Knock-Out of Retrovirus Receptor Gene Tva in the Chicken Confers Resistance to Avian Leukosis Virus Subgroups A and K and Affects Cobalamin (Vitamin B12)-Dependent Level of Methylmalonic Acid.

The chicken Tva cell surface protein, a member of the low-density lipoprotein receptor family, has been identified as an entry receptor for avian leukosis virus of classic subgroup A and newly emerging subgroup K. Because both viruses represent an important concern for the poultry industry, we introduced a frame-shifting deletion into the chicken tva locus with the aim of knocking-out Tva expression and creating a virus-resistant chicken line. The tva knock-out was prepared by CRISPR/Cas9 gene editing in chicken primordial germ cells and orthotopic transplantation of edited cells into the testes of sterilized recipient roosters. The resulting tva -/- chickens tested fully resistant to avian leukosis virus subgroups A and K, both in in vitro and in vivo assays, in contrast to their susceptible tva +/+ and tva +/- siblings. We also found a specific disorder of the cobalamin/vitamin B12 metabolism in the tva knock-out chickens, which is in accordance with the recently recognized physiological function of Tva as a receptor for cobalamin in complex with transcobalamin transporter. Last but not least, we bring a new example of the de novo resistance created by CRISPR/Cas9 editing of pathogen dependence genes in farm animals and, furthermore, a new example of gene editing in chicken.

Repeated MDA5 Gene Loss in Birds: An Evolutionary Perspective.

Two key cytosolic receptors belonging to the retinoic acid-inducible gene I (RIG-I)-like receptor (RLR) family sense the viral RNA-derived danger signals: RIG-I and melanoma differentiation-associated protein 5 (MDA5). Their activation establishes an antiviral state by downstream signaling that ultimately activates interferon-stimulated genes (ISGs). While in rare cases RIG-I gene loss has been detected in mammalian and avian species, most notably in the chicken, MDA5 pseudogenization has only been detected once in mammals. We have screened over a hundred publicly available avian genome sequences and describe an independent disruption of MDA5 in two unrelated avian lineages, the storks (Ciconiiformes) and the rallids (Gruiformes). The results of our RELAX analysis confirmed the absence of negative selection in the MDA5 pseudogene. In contrast to our prediction, we have shown, using multiple dN/dS-based approaches, that the MDA5 loss does not appear to have resulted in any compensatory evolution in the RIG-I gene, which may partially share its ligand-binding specificity. Together, our results indicate that the MDA5 pseudogenization may have important functional effects on immune responsiveness in these two avian clades.

Recombination Marks the Evolutionary Dynamics of a Recently Endogenized Retrovirus.

All vertebrate genomes have been colonized by retroviruses along their evolutionary trajectory. Although endogenous retroviruses (ERVs) can contribute important physiological functions to contemporary hosts, such benefits are attributed to long-term coevolution of ERV and host because germline infections are rare and expansion is slow, and because the host effectively silences them. The genomes of several outbred species including mule deer (Odocoileus hemionus) are currently being colonized by ERVs, which provides an opportunity to study ERV dynamics at a time when few are fixed. We previously established the locus-specific distribution of cervid ERV (CrERV) in populations of mule deer. In this study, we determine the molecular evolutionary processes acting on CrERV at each locus in the context of phylogenetic origin, genome location, and population prevalence. A mule deer genome was de novo assembled from short- and long-insert mate pair reads and CrERV sequence generated at each locus. We report that CrERV composition and diversity have recently measurably increased by horizontal acquisition of a new retrovirus lineage. This new lineage has further expanded CrERV burden and CrERV genomic diversity by activating and recombining with existing CrERV. Resulting interlineage recombinants then endogenize and subsequently expand. CrERV loci are significantly closer to genes than expected if integration were random and gene proximity might explain the recent expansion of one recombinant CrERV lineage. Thus, in mule deer, retroviral colonization is a dynamic period in the molecular evolution of CrERV that also provides a burst of genomic diversity to the host population.

The avian retroviral receptor Tva mediates the uptake of transcobalamin bound vitamin B12 (cobalamin).

The Avian sarcoma and leukosis viruses (ASLVs) are important chicken pathogens. Some of the virus subgroups, including ASLV-A and K, utilize the Tva receptor for cell entrance. Though Tva was identified three decades ago, its physiological function remains unknown. Previously, we have noted an intriguing resemblance and orthology between the chicken gene coding for Tva and the human gene coding for CD320, a receptor involved in cellular uptake of transcobalamin (TC) in complex with vitamin B12/cobalamin (Cbl).Here we show that both the transmembrane and the glycosylphosphatidylinositol (GPI)-anchored form of Tva in the chicken cell line DF-1 promotes the uptake of Cbl with help of expressed and purified chicken TC. The uptake of TC-Cbl complex was monitored using an isotope- or fluorophore-labeled Cbl. We show that (i) TC-Cbl is internalized in chicken cells; and (ii) the uptake is lower in the Tva-knockout cells and higher in Tva-overexpressing cells when compared with wild type chicken cells. The relation between physiological function of Tva and its role in infection was elaborated by showing that infection with ASLV subgroups (targeting Tva) impairs the uptake of TC-Cbl, while this is not the case for cells infected with ASLV-B (not recognized by Tva). In addition, exposure of the cells to a high concentration of TC-Cbl alleviates the infection with Tva-dependent ASLV.IMPORTANCE: We demonstrate that the ASLV receptor Tva participates in the physiological uptake of TC-Cbl, because the viral infection suppresses the uptake of Cbl and vice versa. Our results pave the road for future studies addressing the issues: (i) whether a virus infection can be inhibited by TC-Cbl complexes in vivo; and (ii) whether any human virus employs the human TC-Cbl receptor CD320. In broader terms, our study sheds light on the intricate interplay between physiological roles of cellular receptors and their involvement in virus infection.

Antiviral Activity and Adaptive Evolution of Avian Tetherins.

Tetherin/BST-2 is an antiviral protein that blocks the release of enveloped viral particles by linking them to the membrane of producing cells. At first, BST-2 genes were described only in humans and other mammals. Recent work identified BST-2 orthologs in nonmammalian vertebrates, including birds. Here, we identify the BST-2 sequence in domestic chicken (Gallus gallus) for the first time and demonstrate its activity against avian sarcoma and leukosis virus (ASLV). We generated a BST-2 knockout in chicken cells and showed that BST-2 is a major determinant of an interferon-induced block of ASLV release. Ectopic expression of chicken BST-2 blocks the release of ASLV in chicken cells and of human immunodeficiency virus type 1 (HIV-1) in human cells. Using metabolic labeling and pulse-chase analysis of HIV-1 Gag proteins, we verified that chicken BST-2 blocks the virus at the release stage. Furthermore, we describe BST-2 orthologs in multiple avian species from 12 avian orders. Previously, some of these species were reported to lack BST-2, highlighting the difficulty of identifying sequences of this extremely variable gene. We analyzed BST-2 genes in the avian orders Galliformes and Passeriformes and showed that they evolve under positive selection. This indicates that avian BST-2 is involved in host-virus evolutionary arms races and suggests that BST-2 antagonists exist in some avian viruses. In summary, we show that chicken BST-2 has the potential to act as a restriction factor against ASLV. Characterizing the interaction of avian BST-2 with avian viruses is important in understanding innate antiviral defenses in birds.IMPORTANCE Birds are important hosts of viruses that have the potential to cause zoonotic infections in humans. However, only a few antiviral genes (called viral restriction factors) have been described in birds, mostly because birds lack counterparts of highly studied mammalian restriction factors. Tetherin/BST-2 is a restriction factor, originally described in humans, that blocks the release of newly formed virus particles from infected cells. Recent work identified BST-2 in nonmammalian vertebrate species, including birds. Here, we report the BST-2 sequence in domestic chicken and describe its antiviral activity against a prototypical avian retrovirus, avian sarcoma and leukosis virus (ASLV). We also identify BST-2 genes in multiple avian species and show that they evolve rapidly in birds, which is an important indication of their relevance for antiviral defense. Analysis of avian BST-2 genes will shed light on defense mechanisms against avian viral pathogens.

Insight into the epigenetic landscape of a currently endogenizing gammaretrovirus in mule deer (Odocoileus hemionus).

Endogenous retroviruses (ERVs) constitute a significant part of vertebrate genomes. They originated from past retroviral infections and some of them retain transcriptional activity. The key mechanism avoiding uncontrolled ERV transcription is DNA methylation-mediated epigenetic silencing. Despite numerous studies describing the involvement of ERV activity in cellular processes, epigenetic regulation of ERVs is still poorly understood. We previously described a cervid endogenous retrovirus (CrERV) in the mule deer genome. This virus exhibits massive insertional polymorphism, suggesting recent activity. Here we employed NGS-based strategy to determine the methylation pattern of CrERV integrations in four mule deer. Besides the vast majority of methylated integrations, we identified a tiny fraction of demethylated proviral copies. These copies represent evolutionary older integrations located near gene promoters. In general, our work is a first attempt to characterize the epigenetic landscape of insertionally polymorphic ERV on a whole-genome scale and offers insight into its interactions with a host.

Deltaretroviruses have circulated since at least the Paleogene and infected a broad range of mammalian species.

The Deltaretrovirus genus of retroviruses (family Retroviridae) includes the human T cell leukemia viruses and bovine leukemia virus (BLV). Relatively little is known about the biology and evolution of these viruses, because only a few species have been identified and the genomic 'fossil record' is relatively sparse. Here, we report the discovery of multiple novel endogenous retroviruses (ERVs) derived from ancestral deltaretroviruses. These sequences-two of which contain complete or near complete internal coding regions-reside in genomes of several distinct mammalian orders, including bats, carnivores, cetaceans, and insectivores. We demonstrate that two of these ERVs contain unambiguous homologs of the tax gene, indicating that complex gene regulation has ancient origins within the Deltaretrovirus genus. ERVs demonstrate that the host range of the deltaretrovirus genus is much more extensive than suggested by the relatively small number of exogenous deltaretroviruses described so far, and allow the evolutionary timeline of deltaretrovirus-mammal interaction to be more accurately calibrated.

Avian Expression Patterns and Genomic Mapping Implicate Leptin in Digestion and TNF in Immunity, Suggesting That Their Interacting Adipokine Role Has Been Acquired Only in Mammals

In mammals, leptin and tumor-necrosis factor (TNF) are prominent interacting adipokines mediating appetite control and insulin sensitivity. While TNF pleiotropically functions in immune defense and cell survival, leptin is largely confined to signaling energy stores in adipocytes. Knowledge about the function of avian leptin and TNF is limited and they are absent or lowly expressed in adipose, respectively. Employing radiation-hybrid mapping and FISH-TSA, we mapped TNF and its syntenic genes to chicken chromosome 16 within the major histocompatibility complex (MHC) region. This mapping position suggests that avian TNF has a role in regulating immune response. To test its possible interaction with leptin within the immune system and beyond, we compared the transcription patterns of TNF, leptin and their cognate receptors obtained by meta-analysis of GenBank RNA-seq data. While expression of leptin and its receptor (LEPR) were detected in the brain and digestive tract, TNF and its receptor mRNAs were primarily found in viral-infected and LPS-treated leukocytes. We confirmed leptin expression in the duodenum by immunohistochemistry staining. Altogether, we suggest that whereas leptin and TNF interact as adipokines in mammals, in birds, they have distinct roles. Thus, the interaction between leptin and TNF may be unique to mammals.

The Long Pentraxin PTX3 Is of Major Importance Among Acute Phase Proteins in Chickens.

The expression level of acute phase proteins (APPs) mirrors the health status of an individual. In human medicine, C-reactive protein (CRP), and other members of the pentraxin family are of significant relevance for assessing disease severity and prognosis. In chickens, however, which represent the most common livestock species around the world, no such marker has yet gained general acceptance. The aim of this study was therefore, to characterize chicken pentraxin 3 (chPTX3) and to evaluate its applicability as a general marker for inflammatory conditions. The mammalian and chicken PTX3 proteins were predicted to be similar in sequence, domain organization and polymeric structure. Nevertheless, some characteristics like certain sequence sections, which have varied during the evolution of mammals, and species-specific glycosylation patterns, suggest distinct biological functions. ChPTX3 is constitutively expressed in various tissues but, interestingly, could not be found in splenic tissue samples without stimulation. However, upon treatment with lipopolysaccharide (LPS), PTX3 expression in chicken spleens increased to 95-fold within hours. A search for PTX3 reads in various publicly available RNA-seq data sets of chicken spleen and bursa of Fabricius also showed that PTX3 expression increases within days after experimental infection with viral and bacterial pathogens. An experimental infection with avian pathogenic E.coli and qPCR analysis of spleen samples further established a challenge dose-dependent significant up-regulation of chPTX3 in subclinically infected birds of up to over 150-fold as compared to untreated controls. Our results indicate the potential of chPTX3 as an APP marker to monitor inflammatory conditions in poultry flocks.

Characterization of Chicken Tumor Necrosis Factor-α, a Long Missed Cytokine in Birds.

Tumor necrosis factor-α (TNF-α) is a pleiotropic cytokine playing critical roles in host defense and acute and chronic inflammation. It has been described in fish, amphibians, and mammals but was considered to be absent in the avian genomes. Here, we report on the identification and functional characterization of the avian ortholog. The chicken TNF-α (chTNF-α) is encoded by a highly GC-rich gene, whose product shares with its mammalian counterpart 45% homology in the extracellular part displaying the characteristic TNF homology domain. Orthologs of chTNF-α were identified in the genomes of 12 additional avian species including Palaeognathae and Neognathae, and the synteny of the closely adjacent loci with mammalian TNF-α orthologs was demonstrated in the crow (Corvus cornix) genome. In addition to chTNF-α, we obtained full sequences for homologs of TNF-α receptors 1 and 2 (TNFR1, TNFR2). chTNF-α mRNA is strongly induced by lipopolysaccharide (LPS) stimulation of monocyte derived, splenic and bone marrow macrophages, and significantly upregulated in splenic tissue in response to i.v. LPS treatment. Activation of T-lymphocytes by TCR crosslinking induces chTNF-α expression in CD4+ but not in CD8+ cells. To gain insights into its biological activity, we generated recombinant chTNF-α in eukaryotic and prokaryotic expression systems. Both, the full-length cytokine and the extracellular domain rapidly induced an NFκB-luciferase reporter in stably transfected CEC-32 reporter cells. Collectively, these data provide strong evidence for the existence of a fully functional TNF-α/TNF-α receptor system in birds thus filling a gap in our understanding of the evolution of cytokine systems.

Remnants of an Ancient Deltaretrovirus in the Genomes of Horseshoe Bats (Rhinolophidae).

Endogenous retrovirus (ERV) sequences provide a rich source of information about the long-term interactions between retroviruses and their hosts. However, most ERVs are derived from a subset of retrovirus groups, while ERVs derived from certain other groups remain extremely rare. In particular, only a single ERV sequence has been identified that shows evidence of being related to an ancient Deltaretrovirus, despite the large number of vertebrate genome sequences now available. In this report, we identify a second example of an ERV sequence putatively derived from a past deltaretroviral infection, in the genomes of several species of horseshoe bats (Rhinolophidae). This sequence represents a fragment of viral genome derived from a single integration. The time of the integration was estimated to be 11-19 million years ago. This finding, together with the previously identified endogenous Deltaretrovirus in long-fingered bats (Miniopteridae), suggest a close association of bats with ancient deltaretroviruses.

Male fertility restored by transplanting primordial germ cells into testes: a new way towards efficient transgenesis in chicken.

The ongoing progress in primordial germ cell derivation and cultivation is opening new ways in reproductive biotechnology. This study tested whether functional sperm cells can be matured from genetically manipulated primordial germ cells after transplantation in adult testes and used to restore fertility. We show that spermatogenesis can be restored after mCherry-expressing or GFP-expressing primordial germ cells are transplantated into the testes of sterilized G0 roosters and that mCherry-positive or GFP-positive non-chimeric transgenic G1 offspring can be efficiently produced. Compared with the existing approaches to primordial germ cell replacement, this new technique eliminates the germ line chimerism of G0 roosters and is, therefore, faster, more efficient and requires fewer animals. Furthermore, this is the only animal model, where the fate of primordial germ cells in infertile recipients can be studied.

Characterization of Two Historic Smallpox Specimens from a Czech Museum.

Although smallpox has been known for centuries, the oldest available variola virus strains were isolated in the early 1940s. At that time, large regions of the world were already smallpox-free. Therefore, genetic information of these strains can represent only the very last fraction of a long evolutionary process. Based on the genomes of 48 strains, two clades are differentiated: Clade 1 includes variants of variola major, and clade 2 includes West African and variola minor (Alastrim) strains. Recently, the genome of an almost 400-year-old Lithuanian mummy was determined, which fell basal to all currently sequenced strains of variola virus on phylogenetic trees. Here, we determined two complete variola virus genomes from human tissues kept in a museum in Prague dating back 60 and 160 years, respectively. Moreover, mass spectrometry-based proteomic, chemical, and microscopic examinations were performed. The 60-year-old specimen was most likely an importation from India, a country with endemic smallpox at that time. The genome of the 160-year-old specimen is related to clade 2 West African and variola minor strains. This sequence likely represents a new endemic European variant of variola virus circulating in the midst of the 19th century in Europe.